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■ Abstract 

f*->) . The static three-quark (3Q) potential is studied in SU(3) lattice QCD with 12 3 x 24 and (3 = 5.7 

' at the quenched level. From the 3Q Wilson loop, 3Q ground-state potential V3Q is extracted using 

the smearing technique for ground-state enhancement. With accuracy better than a few %, V3Q 
is well described by a sum of a constant, the two-body Coulomb term and the three-body linear 
confinement term a"3Q_L m i n , with L m i n the minimal value of total length of color flux tubes linking 
the three quarks. Comparing with the Q-Q potential, we find a universal feature of the string 

■ tension, 03Q ~ £>"qq, and the OGE result for Coulomb coefficients, Azq ~ ^Aqq. 

(N ' 

I In usual, the three-body force is regarded as a residual interaction in most fields of physics. In QCD, 

however, the three-body force among three quarks is expected to be a "primary" force reflecting the 
SU(3) C gauge symmetry. Indeed, the three quark (3Q) potential ||, |) is directly responsible to the 
structure and properties of baryons similar to the relevant role of the Q-Q potential upon meson 
properties ||. In contrast with a number of studies on the Q-Q potential using lattice QCD ||, [?[], 
there were only a few lattice QCD studies for the 3Q potential done mainly more than 13 years ago 
H |9|, [h], [Tl| . (In Ref. jll[] , the author only showed a preliminary result on the equilateral-triangle case 
without enough analyses.) In Refs.|| [l(], |TT[ ] , the 3Q potential seemed to be expressed by a sum of two- 
body potentials, which supports the A- type flux tube picture ]i~2| [. On the other hand, Ref.|9j seemed 
to support the Y-type flux-tube picture (2| Q rather than the A-type one. These controversial results 
" may be due to the difficulty of the accurate measurement of the 3Q ground-state potential in lattice 

QCD. For instance, in Refs.|| the authors did not use the smearing for ground-state enhancement, 
and therefore their results may include serious contamination from the excited-state component. 
•iH . The 3Q static potential can be measured with the 3Q Wilson loop, where the 3Q gauge-invariant 

state is generated at t = and is annihilated at t = T, as shown in Fig.l. Here, the three quarks are 
spatially fixed in R 3 for < t < T. The 3Q Wilson loop W$q is defined in a gauge-invariant manner as 

W 3Q = jfabcEa'Vc'Ufufuf (1) 

with Uk = Pexp{ig f r dx^A fl (x)} (k = 1,2,3). Here, P denotes the path-ordered product along the 
path denoted by Tk in Fig.l. Similar to the derivation of the Q-Q potential from the Wilson loop, the 
3Q potential V3Q is obtained as V3Q = — limr^oo ^ ln(H3Q). 

Physically, the true ground state of the 3Q system, which is of interest here, is expected to be 
expressed by the flux tubes instead of the strings, and the 3Q state which is expressed by the three strings 
generally includes many excited-state components such as flux-tube vibrational modes. Of course, if the 
large T limit can be taken, the ground-state potential would be obtained. However, (W^q) decreases 
exponentially with T, and then the practical measurement of (W^q) becomes quite severe for large T in 
lattice QCD simulations. Therefore, for the accurate measurement of the 3Q ground-state potential V3Q, 
the smearing technique for ground-state enhancement [Q, |l3[ |l4| is practically indispensable. However, 
this smearing technique was not applied to the past lattice QCD studies for V 3 q in Rcfs. [|[ |l0), since 
the smearing technique were mainly developed after their works. 

In this paper, we study the 3Q ground-state potential V3Q using the ground-state enhancement by 
the gauge-covariant smearing technique for the link-variable in SU(3) C lattice QCD with the standard 
action with (3=5.7 and 12 3 x 24 at the quenched level jij. We consider 16 patterns of the 3Q configuration 
where the three quarks are put on (i, 0, 0), (0, j, 0) and (0, 0, k) in R 3 with < i, j, k < 3 in the lattice 
unit. Here, the junction point O is set at the origin (0,0,0) in R 3 , although the final result of the 
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ground-state potential V3Q should not depend on the artificial selection of O. For actual lattice QCD 
calculations of the 3Q Wilson loop, we use the translational, the rotational and the reflection symmetries 
on lattices on the choice of the origin O and the direction of x, y, z 

The standard smearing for link- variables is expressed as the iteration of the replacement of the spatial 
link-variable Ui(s) (i = 1,2,3) by the obscured link-variable Oj(s) £ SU(3) C [[l3| ff3 which maximizes 
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aU^s) + ^{Uj&Uiis + jpjis + 1) + Uj(s- j)Ui(s - j)Uj( 8 + I - j)} 
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with the smearing parameter a being a real number. The n-th smeared link- variables E7^(s) (n — 
1, 2, .., iVsmear) are iteratively defined starting from uj^'(s) = Up,(s) as 

utHs) = U^ n - 1 \s) (<= 1,2,3), ut\s) = U i (s). (3) 

As an important feature, this smearing procedure keeps the gauge covariance of the "fat" link- 
variable ulr\s) properly. In fact, the gauge-transformation property of [///^(s) is just the same as that 
of the original link- variable Uu(s), and therefore the gauge invariance of FfJjjT'Hs)) is ensured whenever 
F(Ufj,(s)) is a gauge- invariant function. 

Since the fat link- variable ujT\s) includes a spatial extension, the "line" expressed with ujF\s) 
physically corresponds to a "flux tube" with a spatial extension. Therefore, if a suitable smearing is 
done, the "line" of the fat link-variable is expected to be close to the ground-state flux tube. (This 
smearing technique is actually successful for the extraction of the Q-Q potential in lattice QCD fu(].) 

Now, we investigate the magnitude of the ground-state component in the 3Q-state operator at 
t = 0, T in the 3Q Wilson loop W3Q. From the similar argument in the Q-Q system [jljj, the overlap of 
the 3Q-state operator with the ground state is estimated with 

Co = (W 3Q (T)) T+1 /(W 3Q (T + 1)) T , (4) 

which we call the ground-state overlap. If the 3Q state at t = 0, T in Fig.l is the perfect ground state, 
(W 3 q(T)) = e - v ^ T and then C =l hold. (Here, W 3Q (T) in Eq.(l) is normalized as (W 3Q (T = 0)) = 1.) 

The ground-state potential V3Q can be measured accurately if Co is large enough. Then, we check 
the ground-state overlap Co of the 3Q Wilson loop (W / 3q(C/^™^(s),T)) composed of the fat link- variable 
U < ^ l \s) in lattice QCD simulations, and search reasonable values of the smearing parameter a and the 
iteration number A^ smear for this purpose. In lattice QCD simulations, the ground-state overlap Co is 
largely enhanced as 0.8 < Co < 1 for T < 3 by the smearing with a = 2.3 and N smeal = 12 for all of the 
3Q configurations in consideration, as shown in Fig. 2. 

Here, we make a theoretical consideration on the potential form in the Q-Q and 3Q systems with 
respect to QCD. In the short-distance limit, perturbative QCD is applicable and the Coulomb-type 
potential appears as the one-gluon-exchange (OGE) result. In the long-distance limit at the quenched 
level, the flux-tube picture would be applicable from the argument of the strong-coupling limit of QCD 
|^|, [j], [Fq |, and hence a linear-type confinement potential proportional to the total flux-tube length is 
expected to appear. Of course, it is nontrivial that these simple arguments on the ultraviolet and 
infrared limits of QCD hold for the intermediate region as 0.2 f m < r < 1 fm. Nevertheless, lattice QCD 
results for the Q-Q ground-state potential is well described by 

V Q Q(r) = -^f- + <r QQ r + C QQ (5) 

at the quenched level |l4| . Actually, we measure Vqq from the on-axis Wilson loop with the smearing 
with a = 2.3 and iV smear = 20, and find a good fitting of Eq.(5) with the parameters (^4qq, Cqq, Cqq) 
listed in Table 1. (In general, the adequate values of a and -/V smear depend on the operator.) In fact, 
Vqq is described by a sum of the short-distance OGE result and the long-distance flux-tube result. 

Also for the 3Q ground-state potential V 3 q , we try to apply this simple picture of the short-distance 
OGE result plus the long-distance flux-tube result. Then, the 3Q potential V-jq is expected to take a 
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form of 

V3Q = ~^3Q ^ ] 1 + 0-3Qimin + C 3 Q, (6) 

where £ m in denotes the minimal value of total length of color flux tubes linking the three quarks. 
Denoting three sides of the 3Q triangle by a, b and c, L m in is expressed as 



-(a 2 + b 2 + c 2 ) + — y/(a + b+ c)(-a + 6 + c)(a - & + c)(a + 6 - c) 

_ Z 



(7) 



when all angles of the 3Q triangle do not exceed 2tt/3. In this case, there appears the physical junction 
which connects the three flux tubes originating from the three quarks, and the shape of the 3Q system 
is expressed as a Y-type flux tube |Q, pjj ], where the angle between two flux tubes is found to be 27r/3 

|j. When an angle of the 3Q triangle exceeds 27r/3, one finds L m { n = a + b + c — max(a, b, c). 

Now, we show lattice QCD results for the 3Q ground-state potential. We generate 210 gauge config- 
urations using SU(3) C lattice QCD Monte-Carlo simulation with the standard action with j3 = 5.7 and 
12 3 x 24 at the quenched level. The lattice spacing a ~ 0.19 fm is determined so as to reproduce the 
string tension as ct=0.89 GeV/fm in the Q-Q potential Vqq. Here, the gauge configurations are taken 
every 500 sweeps after a thermalization of 5000 sweeps using the pseudo-heat-bath algorithm. In this 
study, lattice QCD calculations have been performed on NEC-SX4 at Osaka University. 

We measure the 3Q ground-state potential V3Q using the smearing technique, and compare the 
lattice data with the theoretical form of Eq.(6). Owing to the smearing, the ground-state component 
is largely enhanced, and therefore the 3Q Wilson loop (W^q) composed with the smeared link-variable 
exhibits a single-exponential behavior as (W 3 q) ~ e~ VsQT even for a small value of T. Then, for each 
3Q configuration, we extract V3Q * from the least squares fit with the single-exponential form 

(W 3Q ) = Ce~ v ^ T (8) 

in the range of T m ; n < T < T max listed in Table 2. Here, we choose the fit range of T such that 
the stability of the "effective mass" V(T) = \n{(W 3Q (T)) / (W 3 q(T + I))} is observed in the range of 
T ■ < T < T — 1 

-l mm _ J _ ± max x ■ 

For each 3Q configuration, we summarize the lattice data V^q as well as the prefactor C in Eq.(8), 
the fit range of T and xV-^df in Table 2. The statistical error of V^q is estimated with the jackknife 
method. We find again a large ground-state overlap as C > 0.8 for all 3Q configurations. 

Now, we consider the potential form. In Fig. 3, we plot the 3Q ground-state potential V3Q as the 
function of £ mm - Apart from a constant, V3Q is almost proportional to i m ; n in the infrared region. We 
show in Table 1 the best fitting parameters in Eq.(6) for V3Q. On the goodness of this fitting, \ 2 seems 
relatively large as xV-^df = 3.76, which may reflect a systematic error on the finite lattice spacing. 

We add in Table 2 the comparison of the lattice data VjJq* with the fitting function V 3 q in Eq.(6) 
with the parameters listed in Table 1. The three-quark ground-state potential V 3 q is well described by 
Eq.(6) with accuracy better than a few %. 

Next, we compare the coefficients (03Q, ^q, C 3 q) in the 3Q potential V 3 q in Eq.(6) with (ctqq, Aqq, Cq 
in the Q-Q potential Vqq in Eq.(5) as listed in Table 1. As a remarkable fact, we find a universal feature 
of the string tension, <t 3 q ~ ctqq, as well as the OGE result for the Coulomb coefficient, A 3 q ~ 5^4qq- 

As a check, we consider the diquark limit, where two quark locations coincide in the 3Q system. In 
the diquark limit, the static 3Q system becomes equivalent to the Q-Q system, which leads to a physical 
requirement on the relation between V 3 q and V'qq. Our results, 1T3Q ~ Uqq and A 3 q ~ 5^QQi are 
consistent with the physical requirement in the diquark limit. Here, the constant term is to be considered 
carefully in the diquark limit, because there appears a singularity or a divergence from the Coulomb 
term in V3Q in the continuum diquark limit, lim rj _» ri | r _ 3 r Q | = — 00. In the lattice regularization, this 

ultraviolet divergence is regularized to be a finite constant with the lattice spacing a as i^rr^] 



l3Q -a 3 q 



where ui is a dimensionless constant satisfying < uj < 1 and to ~ 1. Then, we find C 3 q H — ^p- = C, 



i.e., C 3 q — Cqq = (> 0) in the diquark limit. This is the requirement for the constant term in 

the diquark limit on the lattice. Our lattice QCD results for C 3 q, Cqq and A3Q are thus consistent 
with this requirement, and one finds u> ~ 0.46. 
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Finally, we also try to fit Kjq with the A-type flux-tube ansatz of V3Q = — Aa J2i<j | r --r-| 
&Aj2i<j \ r i ~ r j\ + Ca, which was suggested in Refs.[|], |l(], However, this fitting seems rather 
worse, because x 2 is unacceptably large as xV-^df = 10.1 even for the best fit with Aa = 0.1410(64), 
a/A = 0.0858(16) and Ca = 0.9344(210) in the lattice unit. As an approximation, V3Q seems described 
by a simple sum of the effective two-body Q-Q potential with a reduced string tension as ctqq ~ 0.53<r. 
This reduction factor can be naturally understood as a geometrical factor rather than the color factor, 
since the ratio between L m ; n and the perimeter length Lp satisfies | < -^j 2 - < -7=, which leads to 
L m in<r = LpUQQ with ctqq = (0.5 — 0.58)<r. 

We have studied the static 3Q ground-state potential in SU(3) lattice QCD at the quenched level, 
using the smearing technique for ground-state enhancement. We have found that V3Q is well described 
by a sum of a constant, the two-body Coulomb term and the three-body linear confinement term 
f3Qimin, where L m i n denotes the minimal length of the color flux tube linking the three quarks. We 
have also found a universal feature of the string tension as 03Q ~ oqq, and the OGE result for Coulomb 
coefficients as A3Q ~ |^4qq- In lattice QCD studies, however, there appear systematic errors relating 
to the finite lattice-spacing effect and so on. To obtain the conclusive result on the 3Q potential, we are 
investigating with finer lattices with large /3's. 
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Figure 1: The 3Q Wilson loop W3Q. The 3Q state is generated at t 
The three quarks are spatially fixed in R 3 for < t < T. 



and is annihilated at t 
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Figure 2: The ground-state overlap of the 3Q system, C = (W 3Q {T)) T+1 / (W 3Q (T + 1)) T , with the 
smeared link- variable (upper data) and with unsmeared link- variable (lower data). To distinguish the 
3Q system, we have taken the horizontal axis as L m i n , which the minimal length of the flux tubes linking 
the three quarks. For each 3Q configuration, Co is largely enhanced as 0.8 < Co < 1 by the smearing. 
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Table 1: The coefficients in Eq.(6) for the 3Q potential and those in Eq.(5) for the Q-Q potential in 
the lattice unit. 





a 


A 


C 


3Q 


0.1524(28) 


0.1331(66) 


0.9182(213) 


QQ 


0.1629(47) 


0.2793(116) 


0.6203(161) 



Table 2: Lattice QCD results for the 3Q potential V^q* for 16 patterns of the 3Q system, where the 
three quarks are put on (i, 0, 0), (0, j, 0) and (0, 0, k) in R 3 in the lattice unit. For each 3Q configuration, 
V^q* is measured from the single-exponential fit (Wsq) — Ce~ V3QT in the range of T listed at the fourth 
column. The statistical errors listed are estimated with the jackknife method, and x 2 /Wdf is listed at 
the fifth column. The fitting function V^q in Eq.(6) with the best fitting parameters in Table 1 is also 
added. 
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Figure 3: The 3Q ground-state potential V3Q as the function of L min which is the minimum value of 
the total length of the flux tubes. 
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